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r~| Abstract 

Gauge mediation provides us with a complete picture of supersymmetry break- 
ing and its mediation within the effective field theories, and thus allows us to 
discuss consistencies with low-energy particle physics as well as cosmological 
Qf) observations. We study in detail the cosmological evolution of the pseudo-modulus 

lO field in the supersymmetry breaking sector and also the production of the graviti- 

nos in the early Universe in a simple (but a complete) model of gauge mediation. 
Under fairly reasonable assumptions, it is found that there exists a parameter 
region where dark matter of the Universe is explained by both thermally and 
non-thermally produced gravitinos, while the baryon asymmetry of the Universe 
. . is generated through the thermal leptogenesis. 
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1 Introduction 



If supersymmetry is realized in nature, there should be a fermionic superpartner of the 
graviton, the gravitino. Since supersymmetry (SUSY) must be spontaneously broken, the 
gravitino is massive. Nevertheless, it becomes massless when we send the Planck scale to 
infinity, and thus the gravitino is the lightest superparticle under an assumption that the 
gravity is the weakest force to talk to the SUSY breaking sector. This discussion already 
motivates us to consider the gravitino as dark matter of the Universe, regardless of the 
hierarchy problem. The assumption that the gravity as the weakest communication is always 
true in gauge mediation scenarios, where various explicit models are known. We take one of 
the simple models of gauge mediation as an example, and demonstrate how such a scenario 
is consistent with observations of the Universe. 

In linearly realized models of SUSY breaking, i.e., in the O'Raifeartaigh models, there is 
a chiral multiplet S whose F-component acquires a vacuum expectation value (VEV). The 
fermionic component of S is the Nambu-Goldstone fermion (the goldstino), and it is eaten by 
the gravitino in supergravity theory. The S multiplet also contains a complex (or two real) 
physical scalar field which radiatively obtains a mass. Since the scalar component has no 
potential at tree level or at the supersymmetric level, it is called the pseudo- modulus, that 
may have cosmological importance. We will consider the dynamics of this field later in detail. 

In the early Universe, the gravitinos are produced in various ways. The scattering of 
particles with superparticles in the thermal plasma, for example, produce gravitinos [TJ El 
El HI El El El El E] • The fraction of the gravitino energy density, n 3 / 2 , from such thermal 
productions is proportional to the reheating temperature of the Universe, TrJ*] Since fi 3 / 2 
is bounded by the total dark matter density ^dm ~ 0.2, Tr needs to be low enough, and 
that possibly causes a tension with scenarios of baryogenesis. For example, the thermal 
leptogenesis [10] requires Tr > 10 9 GeV [TT] which is too high in typical gauge mediation 
models. Gravitinos can also be produced non-thermally such as by the decays of the string 
moduli EH H51 EH 03 HE!, the mflaton [H [201 E0 ES [23] and/or the pseudo-modulus 

field [21 ESI H3 EZ1 Eg El EH]. 

Once we set up a whole picture of a gauge mediation model, the gravitino abundance 
through the thermal scattering processes and also through non-thermal productions via the 
pseudo-moduli decays are calculable. For example, in Ref. [281 [29], a simple model of gauge 
mediation [31] has been considered, and indeed it has been found that the non-thermally 
* Throughout this paper, the reheating refers to the decay of the infiaton. 
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generated gravitino can explain the dark matter abundance in a theoretically motivated 
parameter region. 

In fact, the dynamics of the pseudo- modulus is somewhat complicated. As the theorem 
of Ref. |32j states, the O'Raifeartaigh-type model generally has multiple potential minima, 
and the one which gives the gaugino mass correctly is never the lowest one. Although the 
theorem does not apply to the model of [5Tj . the situation is similar due to the existence of 
a supersymmetric vacuum where the messenger fields condense. An important assumption 
made in the studies of Refs. [28| 29j is that the field value of the pseudo-modulus after inflation 
is located far away from its origin so that it never approaches to the true minimum. Under 
this assumption, it has been found that the pseudo-modulus successfully settles down at the 
meta-stable SUSY breaking vacuum. 

However, it has been pointed out recently that the above assumption is not necessary 
when we take into account the thermal effects on the potential |33| I34j . Once we start 
from the origin of the pseudo-modulus, the messenger fields are massless and thus they are 
thermalized. The thermal effects of the messenger fields contribute to the potential of the 
pseudo-modulus. A parameter region where the meta-stable minimum is selected along the 
thermal history was identified. 

It is then natural to ask whether the gravitino abundance is consistent with observations 
in such a scenario. Since the messenger fields are abundant in the thermal plasma, it looks 
dangerous for the production of the goldstino component of the gravitino which is directly 
coupled to the messenger fields. In order to see if this scenario is viable, one needs to follow 
the cosmological history and calculate the gravitino abundance both from the thermal and 
non-thermal processes. 

In this paper, we explicitly calculate the gravitino abundance in a scenario where the 
initial position of the pseudo-modulus is close to the origin where the thermal potential is 
minimized. We find various non-trivial behavior of the pseudo-modulus depending on model 
parameters. By numerically following the motion of the pseudo- modulus in the field space, 
it is found that the coherent oscillation of the pseudo-modulus eventually dominates the 
energy density of the Universe in a wide range of the parameter space and there the non- 
thermal production of the gravitinos by its decay can explain the right amount of dark matter. 
Moreover, it is found that the reheating temperature after inflation can be much higher 
than the one required by the thermal leptogenesis scenario, without having a trouble with 
the overproduction of the gravitinos. By considering the dilution of the baryon asymmetry 
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by the entropy production from the decays of the pseudo- modulus, the required reheating 
temperature is higher or comparable to 10 12 GeV, with which one can explain both dark 
matter and baryon asymmetry of the Universe. 

The paper is organized as follows. In section 2 we set up a gauge mediation model which 
we use for the study of cosmology. The dynamics of the pseudo-modulus is studied in section 
3 and the results are used in section 4 for the calculation of the gravitino abundance generated 
from scattering processes in the thermal bath. The non-thermal component is calculated in 
section 5. 



2 A model of gauge mediation 



We study the low-energy effective theory of O'Raifeartaigh type SUSY breaking model 
coupled with the messenger fields. After integrating out the massive fields, the Kahler 
potential and the superpotential are written as 

K = Pf + Pf + S^S- {S ^ + --- , (1) 
W = m 2 S - XSff + c, (2) 

where S is a gauge singlet field and / and / are the messenger fields which carry standard 
model quantum numbers. We take the messenger fields / and / to transform as 5 and 5 under 
SU(5), and the messenger number N = 1 for simplicity. A is the cutoff scale of the effective 
theory, which is typically the mass scale of the massive fields of the O'Raifeartaigh models. 
The U(1)r charge assignment is R(S) = 2 and R(ff) = 0, whereas the constant term c 
represents the R-symmetry breaking supergravity effect and contributes to the cosmological 
constant. We take m, A, c and A real and positive by appropriate redefinition of the fields 
and the U(1)r transformation. There is a quantum correction to the Kahler potential from 
the interaction term XSff, 

Kl - loop = ~~^ SiSl ° g A?' (3) 

at the one-loop level. As we shall see, if this radiative correction is too large, the SUSY 
breaking vacuum is destabilized, and so, A is bounded above. 

Once we take into account the supergravity effects, a SUSY breaking vacuum appears [31], 

<S> = ^~, </> = </> = 0, (4) 
6 M p i 
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where we have neglected the radiative correction to the Kahler potential. Here M p \ ~ 2.4 x 
10 18 GeV is the reduced Planck scale. The i^-term of 5* is 

F s = m 2 , (5) 

in this vacuum, and therefore SUSY is indeed broken. The messenger fields acquire a SUSY 
mass M mess = \{S) through the interaction with S. Note that S develops a VEV along the 
real component, while the imaginary component vanishes in this vacuum. The constant term 
c is fixed so that the cosmological constant is cancelled at the SUSY breaking vacuum, 

c ~ m 3/2 Mj, (6) 

where the gravitino mass is given by 

So far we have neglected the radiative correction ([3]), which could destabilize the SUSY 
breaking vacuum. In order to see its effect, let us study the mass terms of the S field and 
the messenger fields (we take / = /), including the correction: 

U mass = (St S)M% (| t ) + (P f)M 2 (p) , (8) 

m 4 / a 15 A 2 M pl \ 

Ml ~ — o A 2 " (9) 

J S ~ A 2 [ 15A 2 ^ 4 j ' 1 ' 

M 2 - ( X2{S)2 ~ Xm2 ) (10) 
M f ~ { -Am 2 A 2 (S> 2 ) ' (10) 

where we have fixed the SUSY breaking vacuum as Q for simplicity, and we have dropped 
terms proportional to A 2 in the diagonal components of M. 2 S - In order for the SUSY breaking 
minimum to be (meta)stable, detmg > and detm 2 > must be met, namely, 



12m 2 M 2 ^ 4^ A 



< A < Ttt^TFT - (11) 



A 4 Vl5 M pl 

Note that A is bounded both above and below. 

There is also a SUSY preserving vacuum at 

(S)SUSY = 0, (/)sUSY = (/)sUSY = \/ (12) 
and the messenger directions are tachyonic for 

\S\<S c ^\f™^. (13) 
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In the absence of the thermal corrections, therefore, S must be outside this region throughout 
the evolution of the Universe since otherwise it would end up with the SUSY minimum. 

Gaugino masses are calculated to be 

A (4tt) 2 (S) (4^) 2 A 2 ' 1 j 

This simple model has three parameters m, A and A. In following analysis we fix the ratio 
m/A so that the gaugino mass should be O(100)GeV. In particular, we fix the Bino mass 
nig ~ 300GeVj^] Once the gaugino mass is fixed, we can relate the gravitino mass to A as 



m 



3/2 



= » 6GeV (3^lv)(l0T4v) 2 - < 15 > 



3 Pseudo-modulus dynamics 
3.1 Vacuum selection 



As the model breaks SUSY at a meta-stable state, we have to check whether the SUSY 
breaking minimum is actually selected in the cosmological evolution. The vacuum selection 
in meta-stable SUSY breaking models has been discussed in several literature for the ISS-type 
models [3S 113 ED H21 S3 Si SS] and the generalized O'Raifeartaigh models [Ml HTJ SHI E3 
134] . In particular, it was shown in Ref. [33] that the pseudo-modulus successfully moves from 
near the origin to the SUSY breaking minimum in the model of Eqs. ([I]) and ([2]), by virtue of 
the finite temperature potentials. Here we follow their discussion and clarify the parameter 
region where S successfully reaches SUSY breaking vacuum. 

We assume that the minimal supersymmetric standard model (MSSM) superfields and 
the messenger superfields are in thermal equilibrium in the early Universe. Although S 
superfield is not in thermal plasma, its scalar potential receives thermal corrections due to 
interactions with the messenger fields. The minimum of the thermal potential is located at 
the origin, (S) = (f) = (/) = 0, rather than the SUSY breaking vacuum in Eq. Thus, 
for a sufficiently high temperature, the potential minimum is close to the origin, and its real 
component will gradually approaches the SUSY breaking one as the temperature decreases, 
whereas the imaginary component remains stabilized at the origin. 

t The SM-like Higgs boson with mass ~ 125 GeV was recently found by the ATLAS \So\ and CMS [HS] 
experiments. For the Bino mass adopted in the text, the light Higgs boson mass may not be explained, unless 
the Higgs sector is extended or extra vector-like matter |371 138] is introduced. Such extensions will not alter 
the following analysis significantly as we mainly focus on the SUSY breaking and mediation sectors. 
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During inflation, on the other hand, there is no thermal plasma. Since the inflaton po- 
tential largely breaks SUSY, the scalar potential of S is modified through Planck-suppressed 
couplings. In particular, the S field generically acquires a so-called Hubble-induced mass, 
and as long as the XJ(1)r remains a good symmetry during inflation, the origin of S is close 
to the extremum of the potential. If the Hubble-induced mass is positive, therefore, S is 
stabilized near the origin during inflation, and it likely remains there even after the inflation, 
which realizes the initial condition mentioned above. 



Note that this initial condition is different from that of Ref. [29], where the S field is 
assumed to be away from the origin in order not to fall into the SUSY preserving vacuum. 
In particular, it was pointed out that the imaginary part of S needs to have a large initial 
value to avoid falling into the wrong vacuum, and the oscillation of the imaginary part is 
eventually the most important because of its relatively long life time. In the present case, 
since the initial position of S is close to the origin after inflation, the imaginary component of 
S stays at the origin throughout, which allows us to concentrate only on the real component 
of S. 

The pseudo-modulus S acquires a finite temperature potential from interactions with the 
messenger fields. The relevant terms along / = / = direction are calculated up to 0(S 3 ): 

V S = -^f~{S + St) + 4^|S| 2 + ^A 2 T 2 |S| 2 - ~A 3 T|S| 3 . (16) 



See Appendix [A] for the derivation. The potential minimum at temperature T is therefore 
given by 



3 A 2 



m 



2 



S min (T)~ V g Mpim 2 + | A2r2 > 



where 



™ S = (18) 

is the tree-level mass of S at zero-temperature. We can see from this formula that the 
potential minimum is near the origin for a high temperature and moves toward the SUSY 
breaking one as the temperature decreases. 

In the absence of the thermal potential, S would fall into the SUSY preserving vacuum if 
\S\ < S a , since the messenger fields become tachyonic. For a sufficiently high temperature, 
however, S does not fall into the SUSY preserving vacuum because the thermal effects lift 
the messenger direction. The messenger fields get thermal potentials through interactions 



S 



with the standard model gauge bosons, 

Vt = -Xm 2 (££ + h.c.) + X 2 \£\ 2 \i\ 2 + X 2 \S\ 2 (\£\ 2 + \£\ 2 ) + ^(3g 2 + g' 2 )(\£\ 2 + \£\ 2 ), (19) 

lb 

V q = -Xm 2 (qq + h.c.) + X 2 \q\ 2 \q\ 2 + X 2 \S\ 2 (\q\ 2 + |g| 2 ) + ^(8g 2 + </ 2 )(|<z| 2 + |g| 2 ), (20) 

where £ and g denote the scalar components of the messenger field /. One can see from the 
potential that the messenger direction becomes unstable at S ~ 0, when the temperature 
becomes lower than T cr : 



T " =4m Vv^' (21) 

As S develops a VEV, the critical temperature goes down. 

We also define a temperature T$ at which S field exits the region where the messenger 
direction is unstable at zero temperature, i.e., 

Smin{Ts) = S CT , (22) 

which gives 

2 8^ m 3 

Ts ~ 15 A3/2M pl ' ( 3) 

In order for the S field to reach the SUSY breaking vacuum without falling into the SUSY 
preserving one, the condition 

T s > T cr (24) 

has to be satisfied^] i.e., the S field should leave the dangerous region (the vicinity of S = 0) 
before the messenger direction becomes tachyonic. This condition is converted to a constraint 
on the model parameters, 

2/5 



A < 



3 3g' z + g' 2 m 



15 16 M 



Pi 



(25) 



We show the allowed region in Fig. [T] From now on, we focus the discussion on the blue 
region where the pseudo-modulus successfully reaches the SUSY breaking vacuum. 



''Precisely speaking, this is a sufficient condition, but we have confirmed that this is consistent with the 
numerical results. 
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13.5 14.0 14.5 15.0 15.5 16.0 16.5 17.0 



Log 10 A[GeV] 

Figure 1: SUSY breaking vacuum (S) appears in the white and blue regions. Starting from 
the origin at high temperature, S field successfully reaches the SUSY breaking vacuum in the 
blue region. The model parameter m is taken to be m = 1.7 x 10~ 7 A so that the Bino mass 
is fixed to be 300GeV. 



3.2 Coherent oscillations 

Let us examine the evolution of S direction more closely. We will find that there are various 
qualitatively different possibilities for the S motion after the inflation. 

During inflation, we assume that S was stabilized near the origin by the positive Hubble- 



induced mass term. After inflation, the S follows the time-dependent minimum (17). When 
the thermal mass becomes comparable to the tree-level mass, the minimum S m m(T) quickly 
moves to the SUSY breaking minimum. This transition takes place at T ~ To, given by 



Whether or not the S field catches up the motion of the minimum depends on the competition 
between the effective mass of S and the friction caused by the expansion of the Universe. 
Here it is assumed that the messenger fields remain thermalized at T = Tq. Later we study 
the case where the messenger fields decouple from the thermal plasma before the temperature 
of the Universe goes down to Tq. 
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The dynamics of S field is governed by the equation of motion 

S + 3HS + -^V = 0. (27) 

The effective mass of S field is approximately given by the sum of the tree-level mass and 
the thermal mass, 

m|(r) = m| + ^A 2 T 2 , (28) 

where we have neglected the contribution from the radiative correction to the Kahler potential 
([3]), since it does not affect results. In the numerical calculations, the effect of the radiative 
correction is properly taken into account. 

First let us consider the case that the Hubble parameter is larger than the effective mass 
at T = To, i.e., H(Tq) > ms(To). In this case, even if the potential minimum moves to the 
zero temperature value at T = Tq, S is still trapped near the origin because of the large 
friction. At a later time when the Hubble parameter becomes comparable to the effective 
mass, H ~ mg, S leaves the vicinity of the origin and starts oscillations about the minimum. 
We define the temperature T osc as 

H{T osc ) = m s {T osc ), (29) 
where the temperature dependence of the Hubble parameter is given by H(T) ~ j^-j- and 



H(T) ~ 2 in the radiation and inflaton-matter dominated eras, respectively. The 



Y n 4 

condition H(Tq) > ms(To) is equivalent to To > T osc . 

On the other hand, if the Hubble parameter is already smaller than the effective mass at 
T = To, H(Tq) < msiTo), or equivalently To < T osc , the friction from the expansion of the 
Universe is small. This is the case if A is larger than the previous case. Then the S field follows 
the time-dependent potential minimum and gradually reaches the SUSY breaking vacuum. 
The amplitude of oscillations is highly suppressed in this case. The suppression was first found 
in Ref . |49| , in which the oscillation amplitude was shown to be exponentially suppressed in a 
limiting cas^] in the context of the cosmological moduli problem. This adiabatic suppression 
mechanism was recently examined more carefully in Ref. [SO]. We show in Fig. [2] the typical 
evolution of S in the above two cases. 

There is yet another possibility. As the temperature of the Universe decreases and the 



value of the pseudo- modulus becomes sizable, the high temperature approximation in Eq.( 16 ) 



§The initial condition adopted in Ref. |49| was given at an infinitely large Hubble parameter. There 
are various additional contributions in general, which are only power-suppressed [50]. We have numerically 
confirmed that the pseudo-modulus abundance is power-suppressed in our scenario. 
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case A 




Figure 2: If ms(To) < H(Tq), S starts coherent oscillations around (S) when H ~ ms (case 
A). The thermal mass increases as A, and the adiabatic suppression takes place (case B). 
The bottom panel show the case of ms(To) > H(Tq), where one can see that S follows the 
minimum S m i n (T) with suppressed oscillations. We set A = 10 16 ' 5 GeV. 

breaks down at a certain point. The messenger fields become non-relativistic when the 
temperature of the Universe becomes comparable to the messenger mass. Then the finite 
temperature potential generated by messenger interactions gets suppressed by the Boltzmann 
factor ~ e~ xs l T . We define the decoupling temperature T^ ec as 

T dec = XS(T = T dcc ). (30) 

If the messenger fields decouple after S field reaches the SUSY breaking vacuum, T^ ec is 
simply the messenger mass scale, 

Tdec ^ A(5) (31) 
for Tdec < Tq. If the decoupling occurs when S is still on the way to the SUSY breaking 
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case C 
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Figure 3: If the messenger fields become non-relativistic before S reaches the zero 
temperature SUSY breaking vacuum (5), the thermal mass term quickly disappears and 
S field begin to oscillate about (S). We set A = 10 165 GeV. 



vacuum, T^ ec is calculated to be 



- dec 



8 a/3 m 4 



1/3 



(32) 



15 AM pl 

for Tdec > Tq. In this case, at T = T& ec , the position of the potential minimum instantly 
moves from S m i n (T) to the SUSY breaking vacuum, which triggers coherent oscillations about 
the minimum (see Fig. [3]). 

In summary, we have defined three temperatures: To, T osc and Td ec - 

• The potential minimum quickly moves from the origin to the SUSY breaking vacuum 
at T = Tq. 

• The Hubble parameter H becomes comparable to the pseudo- modulus mass ms(T) at 
T — T 

± — ^osc- 

• The messenger fields become non-relativistic and disappear from the thermal plasma 
at T = Tci ec . 

The evolution of S, and therefore, its abundance, sensitively depends on the relations 
among these temperatures. The following three cases are shown in Fig.|4]for several reheating 
temperatures. Note that the abundance of S is suppressed in the second case (white region), 
while it is significant in the other cases. 



Tq > T osc (Green region) ; S starts coherent oscillation about (S) at T = T Q: 
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Log 10 A[GeV] 

Figure 4: The evolution of S exhibits distinctive behavior depending on the model 
parameters. In the green region (To > T osc ) S starts coherent oscillations when H ~ ms{T). 
The S abundance is suppressed by the adiabatic suppression mechanism in the white region 
(T osc > Tq > Tdec)- In the blue region (T osc > Td ec > To) the coherent oscillations are 
triggered when messenger fields decouple from thermal plasma and disappear. While To and 
Ti ec are uniquely determined once we fix the model parameters, the value of T osc depends on 
Tr. 



• T osc > Td > Tdoc (White region) ; S follows S m i n (T) and gradually reaches (S) without 
sizable oscillations. 

• T osc > Td cc > Tq (Blue region) ; The messenger fields decouple from thermal plasma 
when S is on the way to (S). Coherent oscillations are triggered at T = Td ec - 
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4 Thermal production of the gravitinos 



In the framework of gauge mediation the dominant component of thermally produced grav- 
itino is that of the longitudinal mode, the goldstino. The relic abundance is estimated to 

be [DUEIEIEIEIE] 



l 3/2, 

where Tr represents the reheating temperature after inflation, and it is assumed that there 



is no entropy production after the reheating. Note however that it is assumed implicitly in 
deriving this formula that the messenger mass scale is higher than Tr. For Tr higher than 
the messenger scale, one needs to consider diagrams in which the messenger fields are in 
the external lines. As shown in Ref. [12], if Tr is higher than the messenger mass scale, or 
equivalently if the messenger fields once gets thermalized, the goldstino relic abundance is 
determined by the messenger mass scale M mess = \{S) rather than the reheating temperature, 

The replacement of Tr by M mess can be understood by looking at the temperature dependence 
of the goldstino reaction rate T(T): 

r(T)~a 3 A 2 T, for T > M mess (35) 

2 

m~ „ 

T(T)~a 3 2 9 2 T for T < M mess . (36) 

m 3/2 M pl 

As r(T) has lower power dependence on T than the Hubble parameter (H oc T 2 ) for 
T > Af mess , contributions to the goldstino production from MSSM fields decouple for T > 
M mess [j] We should estimate the thermally produced gravitino abundance by Eq. ( 34 ) since 



the pseudo-modulus was at the origin and the messenger fields were inevitably thermalized 
in our scenario. 



^ It has been argued in Ref. [T3] that there is a component of the interaction rate which still grows as 
T 3 at high temperatures, and thus the estimation of the gravitino abundance is sensitive to Tr in contrast 
to the conclusion of Ref. 12 . This is based on the observation that there are contact interactions between 
the gravitino and visible sector fields in the supergravity Lagrangian in the unitary gauge. However, we 
anticipate a cancellation of such contributions with loop diagrams at a high temperature since there is no such 
contribution in the goldstino picture. 
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5 Pseudo-modulus decay and the gravitino abundance 



In this section we examine the decays of the pseudo- modulus. So far we have assumed that 
there is no additional entropy production after reheating. If the oscillation energy of the 
pseudo-modulus dominates the Universe, however, a large amount of entropy is produced by 
the decay, and the pre-existing gravitino is diluted. Also the gravitinos are produced by the 
decay of the pseudo- modulus. 



5.1 S'-domination 



The energy density of S decreases more slowly than radiation, and so, if its lifetime is 
sufficiently long or the reheating temperature is high, S may come to dominate the Universe 
and produce a large amount of entropy by the decay. Let us define the entropy dilution factor 
A as 

1 «inf 



- Min 1,^ , (37) 
A S S + Sinf [_ S S J 

where s in f and ss represent the entropy density produced from the inflaton and the pseudo- 
modulus, respectively. If A > 1, the pre-existing gravitino abundance is diluted by a factor 
of 1/A. In this case, we can express the dilution factor in terms of the decay temperature of 
S, T^, and its abundance as 

s s 4 p s 



•Sinf 3Td Si n f 



(38) 

S decay 



The pseudo-modulus abundance, (ps/smi), can be estimated by following its evolution as 
discussed in Sec. |3j The decay temperature Td is defined as 



90 \ 1/4 

T d = [ — ) \/M pl T s , (39) 



7T ff* 



where is the decay width of S. 



The decay width of S can be explicitly calculated. The S field couples to the MSSM 
particles through loop diagrams of the messenger fields. The partial decay width of the two 
gravitino mode is suppressed compared to the MSSM particles. The interaction Lagrangians 
between S and MSSM fields and calculations of decay widths can be found in Refs. [2H| 129"]. 
The main decay mode is S — > bb and S — > hh + WW + ZZ for 2m;, < mj < 2mw and 2m w < 
ms < ITeV cases, respectively. The decay temperature is calculated to be 0(0.1 — 10)GeV 
in the parameter region of our interest. (See Refs. |28t |2"9"] for formulae and parameter 
dependences.) 
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Figure 5: The contour plot of the dilution factor A, for Tr = 10 5 , 10 7 , and 10 12 GeV. The 
colored regions correspond to the regions where the pseudo-modulus dominates the Universe 
before the decay. In the case of Tr = 10 12 GeV, the oscillations of S commence after the 
reheating, and the S tends to dominate the Universe, producing larger entropy compared to 
the other two cases with lower Tr. 



We evaluate the dilution factor numerically. We solve the equation of motion of the 
pseudo- modulus which provides the pseudo- modulus abundance, ps/s- m {. For the potential 



at finite temperatures, we use the one in Eq. (44). The decay temperature is calculated by 
using the formulae in Refs. |28| 129] . The results are shown in Fig. [5] The colored regions 
(A > 1) correspond to the regions where the oscillation energy dominates the energy density 
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of the Universe. One can see that the S abundance is indeed suppressed in the middle of 
each panel, because of the adiabatic suppression mechanism. The similarity with Fig. [4] is 
clear especially in the case of Tr = 10 12 GeV. 

5.2 Non-thermally produced gravitino 

Gravitinos are produced non-thermally by the rare decay S — > ^3/2^3/2 • The branching ratio 
of the gravitino mode, -B3/2, is calculated in Refs. [281 [29], an d it is typically of 0(1O -10 — 10~ 6 ) 
in the parameter regions of our interest. 

If S dominates the energy density of the Universe, non-thermal gravitino abundance is 
calculated as 

^3/ T 2 = ^3/2 — ^3/2 X 2/( Pc / S ) , (40) 

where (p c /s)o ~ 3.6 x 10~ 9 /i~ 2 GeV. We will see later that there are parameter regions where 
the DM relic density is explained by the non-thermally produced gravitinos. 

5.3 Total gravitino abundance 

The relic density of the gravitinos is given by the sum of the thermally and non-thermally 
produced gravitinos. If the S dominates the Universe before the decay, we need to take 
account of the entropy factor, which dilutes the pre-existing gravitinos produced at the 
reheating; 

^3/2 = ^^3/2 + ^3/2- ( 41 ) 

Here the dilution factors are evaluated numerically as in Fig. [5} Note that is so low that 
the thermal production of gravitinos is negligible at the pseudo-modulus decay. 

We show in Fig. [6] the contours of the total gravitino abundance for Tr = 10 5 , 10 7 
and 10 12 GeV. Considering 0(1) uncertainties in the calculations of dilution factor and 
thermal gravitino productions, we show the highlighted region (in red) where the gravitino 
abundance 0.1 < < 1 is obtained. We expect that somewhere in this region will provide 
the observationally consistent DM abundance, I^dm = 0.2. For Tr = 10 5 GeV, the relic 
gravitino is mostly thermally produced one. On the other hand, for Tr = 10 12 GeV, there is 
an allowed region where the dilution factor A is about 10 3 — 10 6 . In this region, the observed 
DM abundance is a mixture of the thermally and non-thermally produced gravitinos. It 
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is crucial to take account of the fact that the gravitino abundance is independent of the 



reheating temperature as in Eq. (34) 



We emphasize here that the contours of the gravitino abundance remain almost intact for 
Tr > 10 9 GeV. This is because the gravitino production rate is modified at a temperature 
above the messenger scale. This has a crucial impact on the leptogenesis as we shall see 
below. 



The gravitino mass in the shaded region ranges from 10 MeV to a few GeV (See Eq. ( 15 )). 
The thermally produced gravitinos therefore behave as cold dark matter. On the other hand, 
the non-thermally produced gravitinos have a relatively large velocity, which results in a 
free streaming length of O(10 — 100) kpc depending on the parameters |28} I29j . In our 
scenario, the DM consists of a mixture of thermally and non-thermally produced gravitinos 
for Tr > 10 7 GeV, and such partial suppression of the density perturbation at small scales 
may have an interesting impact on the large scale structure of the Universe. 

5.4 Leptogenesis 

Interestingly, in this model, it is possible to create the right amount of the baryon asymmetry 
through the thermal leptogenesis |10j while satisfying < 0.2. 

In the thermal leptogenesis, there is an upper bound on the baryon asymmetry for a fixed 
reheating temperature [51j [52j [11] : 

»S tx ^ a < 42 > 

This often causes a tension with the gravitino overproduction. In the presence of the late- 
time entropy production, one can increase the reheating temperature because the thermally 
produced gravitinos are diluted . However, this does not solve the tension because both 



the baryon asymmetry and the conventional formula of the graviton abundance (33) are 
proportional to Tr. 

In our scenario, the reheating temperature can be much higher than 10 9 GeV without 
having a problem of gravitino overproduction, because the gravitino abundance becomes 
independent of the reheating temperature for sufficiently high Tr. Recall that the contours 
of the gravitino abundance in Fig. 6 are almost the same for Tr > 10 9 GeV. 



By comparing Fig. [5] with Fig. [6j we see that the size of the dilution factor A is about 
~ 10 3 — 10 6 in the region where ^3/2 — 0.1 — 1. Therefore, if the reheating temperature is 
higher than ~ 10 12 GeV, even though it is partially diluted, the observed amount of baryon 
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T R = 10 5 GeV T R = 10 7 GeV 
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Log 10 A[GeV] Log 10 A[GeV] 

T R = 10 12 GeV 




14.0 14.5 15.0 15.5 16.0 16.5 17.0 

Log 10 A[GeV] 



Figure 6: The contour plot of the total gravitino abundance Red regions correspond 

to DM consistent regions 0.1 < ^3/2 < 1. In T R = 10 5 GeV case thermally produced gravitino 
comes dominantly from the MSSM sector. The contribution from the messenger sector 
becomes comparable to that of the MSSM sector in T R = 10 7 GeV. In T R = 10 12 GeV case the 
thermal component consists of the messenger sector contribution. Non-thermal components 
from the decay of the pseudo-modulus field are also included in the figures. 

asymmetry can be generated by thermal leptogenesis. This is one of the distinctive features 
of our scenario. 
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6 Summary 



We have followed the cosmological evolutions of the SUSY breaking pseudo-modulus field 
and the messenger fields in a simple gauge mediation model which break SUSY at a meta- 
stable state. We adopt an initial condition such that the pseudo-modulus field is close to the 
origin, stabilized by thermal potential via interactions with the messenger fields in thermal 
plasma. Such an initial condition can be naturally realized if the U(l)# symmetry remains 
a good symmetry during inflation and S is stabilized near the origin due to the positive 
Hubble-induced mass. We have found that this simple gauge mediation model with the 
initial condition is cosmologically viable in a sense that the pseudo-modulus can settle down 
at the correct SUSY breaking minimum and that the gravitino relic abundance can explain 
the DM of the Universe. Furthermore, thermal leptogenesis is possible without the gravitino 
overproduction. 

We have numerically calculated the total gravitino relic abundance, both the thermally 
produced one taking account of the entropy dilution and the non-thermally produced one, 
and shown that the observed DM density can be explained for m 3 / 2 = 0(lO)MeV — C(l)GeV. 

The messenger fields play a crucial role in the scenario. The messengers acquire a thermal 
mass when the pseudo-modulus stays near the origin, which prevents the messengers to fall 
into the SUSY preserving minimum. At the reheating, the gravitino production rate is 
modified if the messenger fields are in the thermal plasma. This has a great impact on the 
gravitino abundance; it becomes independent of the reheating temperature of the Universe. It 
is this fact that enables thermal leptogenesis to create the right amount of baryon asymmetry 
without overproduction of gravitinos. 
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A Pseudo-modulus potential at a finite temperature 



The finite temperature effective potential up to one-loop is given bvji53| 154"] 

V = Vtree + Vi + thermal, (43) 

where Vt iee is the classical potential calculated from Eq.Q and ^ and V\ is the zero- 
temperature one- loop potential calculated from Eq.Q. Finite temperature one- loop correc- 
tion is 



^thermal ^ 2 



/ fa 2 Vio g [i- e -v^K^iy^] 

Jo , 



fOO . POO . 

-2 / dxx 2 V log[l + e-V^+WDr/T 2 ] + 3 / dxx 2 V log[l - e~ ^ 
Jo Jo 



2+(M*) a /Th 



(44) 



where the three terms represent the contributions from real scalar fields fa, Weyl fermions ip r 
and vector bosons A% with the eigenvalues of the squared mass matrices (M|)j, (M|) r and 



(My-) a . For a high temperature limit of T 3> Mg, Mp, My, the potential can be expanded as 



thermal = (n b + ^f) + ^ [Tr(Mf ) + 3Tt(M 2 ) + Tr(M|)] 

-^[!V(Mj) + 3Tr(il4 )]+-... (45) 
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